Asymptotically optimal quantum channel reversal for 
qudit ensembles and multimode Gaussian states 



Peter Bowles, Madalin Guja and Gerardo Adesso 

School of Mathematical Sciences, The University of Nottingham, University Park, NG7 2RD 
Nottingham, United Kingdom 



PACS numbers: 03.67.Hk, 03.65.Wj, 02.50.Tt 



Abstract. We investigate the problem of optimally reversing the action of an arbitrary 
quantum channel C which acts independently on each component of an ensemble of n 
identically prepared rf-dimensional quantum systems. In the limit of large ensembles, we 
construct the optimal reversing channel R* which has to be applied at the output ensemble 
state, to retrieve a smaller ensemble of m systems prepared in the input state, with the highest 
possible rate m/n. The solution is found by mapping the problem into the optimal reversal 
of Gaussian channels on quantum-classical continuous variable systems, which is here solved 
as well. Our general results can be readily applied to improve the implementation of robust 
long-distance quantum communication. As an example, we investigate the optimal reversal 
rate of phase flip channels acting on a multi-qubit register. 
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1. Introduction 

Quantum channels are completely positive, trace preserving maps which describe the state 
change of a quantum system undergoing a noisy evolution. Operationally, a channel on a d- 
dimensional system C : M(C d ) — > M(C d ) can be implemented by first coupling the system 
with state p to the "environment" with initial state p E , letting the two evolve together as 
a closed system with unitary operator U, and subsequently tracing out the environmental 
degrees of freedom (HE) 

C:p^Tr E {U(p®p E )U*). (1) 

Due to the last step, the dynamics described by a quantum channel is typically irreversible, a 
basic fact which is related to several no-go theorems in quantum information [3] [4] [5] |6) . In 
fact, a channel is reversible on all states if and only if it is unitary i.e. C(p) = UpU*. 

However, the reversibility problem has non-trivial solutions if the channel is required to be 
reversible on a given family of states From a mathematical perspective, this scenario is 
captured by the the statistical concepts of quantum sufficiency and equivalence of statistical 
models. In the latter, two families of states 

Q := \p g : 9 e 0} and K := {cr g : e 0} 

are said to be statistically equivalent if there exist two channels T and S such that erg = T(pg) 
and pe — S (erg) for all 6. In other words the channel T is reversible on Q, with an inverse 
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Figure 1. (Color online) Graphical depiction of the problem addressed in this paper, (a) 
Reversal of arbitrary channels C acting independently on a register of n qudits, each initialized 
in the state p. The optimal reversing channel R n allows to recover m qudits exactly in the 
state p, with the maximum possible output/input rate m/n. (b) The problem can be recast into 
the reversal of a Gaussian channel G acting on a set of Gaussian states with unknown mean 
M. The optimal reversing channel H restores the initial known covariance and maximizes the 
output/input displacement ratio k. More details are provided in the diagram |5J and in the 
accompanying text. 



S. Sufficiency is a special case of this, where Tip) is the restriction of the state p to a 
subalgebra of observables, which is then called sufficient with respect to <2. A general 
mathematical characterization of equivalence and sufficiency has been given in Q, where 
a quantum analogue of the classical factorization theorem for sufficient statistics has been 
established. 

In this paper we consider a different but related channel reversal problem, see Fig.[T] Consider 
an ensemble consisting of n » 1 independent and identically prepared quantum systems 
which undergo separate noisy evolutions described by the channel C. The question is whether 
the original state can be distilled from the noisy output by means of a reversing channel, 
albeit at the expense of reducing the size of the ensemble due to the loss of information. More 
precisely, we would like to find the maximum rate T = m/n for which there exist reversing 
channels R„ : M(C d ) m 1— » M{C d )® m such that the following holds 

p m t— > C(pf" hA p ® m (2) 

asymptotically with n (see section|4]for the mathematical formulation). This and other related 
questions have been considered in [4 8 1 for the specific case of a qubit depolarizing channel, 
but to our knowledge the case of a general qudit channel has not been investigated elsewhere. 
The answer to our question is directly relevant to applications in quantum memories ll9l [T0l . 
distributed quantum computation ATI , quantum key distribution [12|, and long-distance 
quantum communication lfl3l . Given a register of many qudits prepared off-line, in order 
to securely store it or transmit it over lossy media, one needs to incorporate modules that 
perform state purification or counteract decoherence effects, e.g. quantum repeaters [ 14 1. 

The key to solving the problem lies in the quantum extension of a fundamental statistical 
tool known as local asymptotic normality (LAN) lfT31 [161 . The quantum version of LAN 
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ifTTl [T8l [T9l l20l Ell shows that the joint state of an ensemble of identically prepared systems 
can be approximated by a Gaussian state of a continuous variable system. More precisely, 
we parametrize the states in a neighbourhood of size n~'/ 2+e of p as Pm/V« wnere M is a 
(d 2 - l)-dimensional vector of expectations such that po = p; then there exists a channel T„ 
which maps the collective state p^ ^_ into a Gaussian state €>(M, V) of a quantum-classical 
continuous variable system, with mean M and known fixed covariance V which depends only 
on p. Conversely, given the Gaussian state, there exists a channel S„ which converts it into the 
multiple-qudits state with asymptotically vanishing norm-one error, without the knowledge of 
the parameter M. 

The schematic of the complete problem considered in this paper and its resolution using LAN 
is shown below (see Fig.[T]for a visual summary): 

m <8« ®m 

^(M,A) ^(M'.A') r (M ,A) 

(3) 



■I -I 



<X>(M,V) — — > 4>(M',V) — ^ 0>(/tM,V) 

Using the LAN correspondence between the multiple-qudits states and the Gaussian states, 
we demonstrate specifically the following results: 

• The action of C m on the qudits can be mimicked by the action of a Gaussian channel on 
the Gaussian states, i.e. the rightmost square loop of the diagram Q is (asymptotically) 
commutative (cf. Theorem |4.1| i; 

• The optimal qudit reversal problem for the channel C is effectively recast into a much 
handier Gaussian one, of reversal of the Gaussian channel G fl22l l23l l24l l25l l26l where 
the length of the final displacement vector is allowed to be a fraction k of the initial length 
as illustrated in the bottom row of the diagram Q; the problem is to find the maximum 



possible k and the corresponding reversing channel H (cf. Theorem 4.2 1; 



The Gaussian reversal problem has an explicit solution (k* ,H*) as described in Theorem 



2.1 The optimal reversal of a general Gaussian channel is in its own right an important 
problem for continuous variable quantum communication and key distribution Il2"8ll30l . 
which we solve in a specific case of input states associated with qudit ensembles. 

• Putting together the above findings we prove that the concatenation R„ := S m o H o T„ 
realizes the optimal strategy for reversal of the channel C m applied to n independent and 
identically prepared qudits, and that the maximal qudit rate m/n in the top row of the 
diagram ^ is equal to (k*) 2 , where k* is the maximum possible value of k for which the 
Gaussian reversal) is possible. 

The paper is organized as follows. In section [2] we find the optimal procedure and rate for 
reversing the action of general Gaussian channels on a family of Gaussian states with fixed 
covariance and unknown mean. Section [3] contains a short review of LAN motivated by the 
quantum Central Limit Theorem. In section [4] we solve the problem of finding the optimal 
rate and procedure for the reversal of a qudit channel by reducing it to the Gaussian problem 
via LAN. In section [5] we illustrate the general results with a specific application to reversing 
phase flip channels acting on a multiqubit register, a common source of errors arising in 
quantum computation HI. We draw our concluding remarks in section [6] Some technical 
proofs are deferred to Appendixes. 
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2. Optimal reversal of a Gaussian channel 

In order to solve the qudit channel reversal problem, it is necessary to obtain the solution to 
the corresponding Gaussian problem. We begin with a short review of Gaussian states and 
channels, referring the reader to l27l [28l [29l l30ll for more comprehensive accounts. 

2.1. Gaussian states and channels 

Let Z = (A, B) r be the (column vector of) coordinates of a continuous variable system where 
the components A r = (Qi,P\, ■ ■ ■ , Qk,Pk) are the canonical observables of K modes, and 
B r = (B\, . . . ,B C ) is a classical R c -valued random variable. In other words, the canonical 
variables satisfy the commutation relations [Zt,Z/] = zO^/l where Q. is the block-diagonal 
symplectic form 

Q. - Diag (cr, . . . ,cr, Oc) , with cr = 

Throughout, we consider that the quantum variables are represented in the standard fashion 
on the multi-mode Hilbert space 

<H®K ii^fK and the c i assica i variables are realized as 
coordinate multiplication operators on < 7/® c so that the full Hilbert space is <J-(®( K + C \ Any state 
<t> of the continuous variable system is completely determined by its characteristic function 

where <&(X) denotes the expectation of X with respect to <I>. In particular, for every mean 
M = 3>(Z) and covariance matrix V with elements 

V u = ®(Z k * Z{) - <fc(Z,)0(Z,), Z k * Z, : = Z * Z/ + Z ' Z * , 

there is a unique Gaussian state €>(M, V) with characteristic function 

XMy (z) = e^-^\ (4) 

provided that V satisfies the uncertainty principle V > - ^Q. For later use, we will denote by 
0(M, V) € T l ( ( H K )) ® L\R C ) the density matrix of the state <D(M, V). 

A Gaussian quantum channel is a channel as in ([TJ, where the environment is a bosonic 
continuous variable system whose initial state is Gaussian, and the unitary U is determined 
by a quadratic Hamiltonian in the system and environment coordinates [22 23 24l l25ll26l . 
In the Heisenberg picture, the action of a Gaussian channel G is 

G* : £(<H® K ) ® L 00 (R l8C ) -> SCH 8 *) ® L°°(R® C ) 

G*:W Z ^ Wx at e-& Yt * (5) 

where W z := e' zTz are the "Weyl operators", Xc, Yq are real matrices of dimension 2K + C, 
with Yc positive and satisfying the matrix inequality 

Y G > i(X T c QX G - Q). 

In particular, from |4]) and |5]l we find that when G acts on the state <t>(M, V), it produces a 
new Gaussian state fl>(M',V) with mean M' = X^M and covariance V = X^VX G + Y c . 
While the first term in V describes the change in variance due to the linear transformation 
Xq, the second term comes from the covariance of the ancillary "environment" used to realize 
G. For more details on Gaussian channels and their classification we refer to 13TI . 
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2.2. Optimal Gaussian channel reversal 

Our problem can be stated as follows: given a Gaussian channel G, and a family of Gaussian 
states 

Qv := {<t>(M, V) : M e R 2 ^ +c ) (6) 

with fixed covariance matrix V and unknown mean M, we would like to find the maximum 
value k* of a constant k for which there exists a (not necessarily Gaussian) channel H such 
that for all M e R 2 ^ +c the following holds [see Fig.[l}b)] 

H oG(0(M,V)) = 4>(kM,Y). (7) 

A special case of this problem, for single-mode attenuation and amplification channels, has 
been considered in [8|. Our first main result is to find k* and the optimal channel for the 
general set-up. 

Theorem 2.1. Let G be the Gaussian channel Q, and assume that Xq is an invertible matrix. 
The largest value of k for which G can be reversed on the family of states (|6} is 

k* = [^(fW^jf' (8) 
where /i max (A) denotes the maximum eigenvalue of A., and 

q = (X G ') r (x£VX G + Y G + -njXe 1 , s = V + 
are positive matrices. 

Equivalently, k* can be expressed in terms of the 'max-relative entropy' H32\l as 

log 2 (fc*~ 2 ) = D m , x (g\W) ■= log 2 (min{r 2 : q < r 2 <r}) . 
For every k < k* there exists a reversing Gaussian channel H, of the form |5]), with 



Y H = V - k 2 (X c ') T (X T c VX c + Y C )X G 1 . 



(9) 



Proof. Following a standard argument [38, 8|, it can be shown that without loss of generality 
we can restrict our attention to a certain class of displacement covariant reversing channels of 
the form 

T*(W Z ) = f(z)W Xz , 

where X = kX c l and f(z) - Tt(tW z ) is the characteristic function of some ancillary state t 
with zero mean. A displacement covariant channel H satisfies (|7| if and only if 

H oG(0(O,V)) = 0(0, V). 

In terms of the characteristic functions, this means 

= ti(H(e izTz )4>') = Tr(TW,)Tr(0y' (X " z)Tz ) 

_ ^ z ) e -i(X«) I- V'(X«z)_ 
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Since x<s> is a Gaussian characteristic functions, f(z) = tr(rW~) must also be a Gaussian 
characteristic function. Hence H is a Gaussian channel. 

Now, from <|3j we see Gaussian channels are in one-to-one correspondence with two matrices: 
the linear transformation on the means and the covariance (noise) matrix. The Gaussian 
channel G and the reversal H are therefore completely characterized by G o (Xq, Yg) and 
H o (Xh, Y//). For an initial state $(M, V), we have that the first and second moments are 
mapped by (H o G) as 

(M, V) ^ (X£X C M, X T H (X T C \X C + Y C )X H + Y H ) 

By equating this with the target state moments (AM, V), we obtain immediately Xh = kX c l . 
Now, denoting V' = X^VXc + Yg, we have that the reverse outputs a covariance 

V = X T H Y'X H + Yh — A^Xc'fV'Xc 1 + Y H . (10) 

Positivity of the quantum channel states that 



Y H +-Q--A 2 (X G 1 ) r QX c 1 >0 (11) 
We would like to find the maximum A for which there exists some Y H which satisfies this 



inequality. Rearranging ( 10 1 as 

Y H =Y-k 2 (X c l ) T YX G [ 

then substituting in ( 1 1 1 gives 

V + > k\X c l ) T (Y + l -Q.)X c l (12) 

This gives us a necessary and sufficient condition for a Gaussian channel G to be reversible 
up to a factor A in the displacement. The optimal reversing channel H* is the one for which A 



takes its maximum possible value A*. To find A*, we can recast ( 12 1 in terms of a max-relative 
entropy and use the results of [32] to obtain Eq. ([8]). □ 

This result shows that any Gaussian channel G acting on Gaussian states with given 
covariance, can be reversed up to a constant factor A by means of another Gaussian channel 
H, whose construction is provided explicitly. The link between the threshold value A* and 
the max-relative entropy l32l is intriguing, as it reveals that A* is operationally related to 
the optimal Bayesian error probability in determining which covariance matrix our system is 
mapped into by the channel G [33 1. 



3. Quick review of Local Asymptotic Normality 

The main tool used to solve the channel reversal problem on n-qudit ensembles in the 
next section of this paper is that of quantum local asymptotic normality (LAN) which was 
developed in ifTTl [181 [T9l l20l I2TI as an extension of a key concept from classical asymptotic 
statistics |[T5l . For reader's convenience we give here an intuitive explanation of LAN based 
on the quantum Central Limit Theorem. 

In general terms, classical LAN means that given data consisting of n samples from a 
probability distribution Pg with unknown parameter 6 € R*, there exist classical channels 
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(randomizations) which map the data into a single sample from a Gaussian distribution 
whose mean is (locally) equal to 9, and whose variance is the inverse Fisher information 
matrix /(#) 1 ■ A consequence of this is the fact that the maximum likelihood estimator 
6„ is asymptotically normal (Gaussian) with asymptotic variance I(6)~ l , i.e. the following 
convergence in distribution holds 



V^(#„ - 0) 



N(0,I(6y l ) 



and therefore saturates the Cramer-Rao bound [16]. 

Quantum LAN means that the joint quantum state of independent, identically prepared finite- 
dimensional systems can be approximated in a strong sense by a quantum-classical Gaussian 
state of fixed variance, whose mean encodes the information about the parameters of the 
original state. In this way, a number of asymptotic problems can be reformulated in terms 
of Gaussian states, for which the explicit solution can be found. Examples so far include 
state estimation [39], teleportation benchmarks [37], state purification and dilution [8| and 
quantum learning lT34ll . 

Suppose we are given n independent li-dimensional quantum systems (qudits) each prepared 
in the unknown but mixed (full rank) state p e M(C d ) with distinct eigenvalues. By means of 
an adaptive estimation strategy, we can effectively localize the initial state of each qudit within 
a neighborhood of size n~~- +e centered at an estimate p ( >. If Ai > • • • > A^ are the eigevalues 
of po then with respect to its eigenbasis, any neighbouring state can be written as 



P(M/ v»U) ~ PO ' 



\fn 



Al,2Zl,2 "2 
AwZl.d 



(13) 



where M = (z, u) e c d(d ~ l)/2 x R d ~ l is a displacement parameter, and A Jk = y/(Aj - A k )/2 
are constant coefficients chosen for later convenience. We can then define the local quantum 
statistical model around po as 



a, 



^) :||M|1 ^' 



Let 



{Ol,....,O d l-\} = fa?, Pi?, -;qd-l,d, Pd-l,d,b U -;b d -l} 

be the selfadjoint operators in M(C d ) defined as 



qj,k := 



\j)(k\ + \j)(k\ 



Pj,k ■= 



i(\k)(j\ - \}){k\) 



^/2(Aj - A k ) ' rJ "' y/2(Aj-A k ) ' 
One can verify easily that O a satisfy the following properties: 



bi-.= \i)(i\-Ail. 



(14) 



(15) 



(16) 



(i) 0j2_i} is a basis in the space of operators with Tr(p o 0) = 0; 

(ii) Tr(po[0fl, Ob]) = iQ a ,b where Q is the (d 2 - l)x(d 2 - 1) block diagonal symplectic matrix 
Q = Diag(cr, . . . , cr, 0^-i). 

(iii) The covariance matrix V with 



V ah := Tr(p o fl * O b ) 



(17) 
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has all elements equal to zero except 
Tv(p q%) = Tv( PoP 2 jk ) = vj* := 2{ i + _l k y ^(Pobjh) := := ^ - AjA k . (18) 
Let (9,(n) 6 M(C')®" denote the corresponding collective (fluctuation) observables 

<9,(n) := Of\ Of := 1 ® • ■ • ® 0,- ® • • ■ ® 1, (19) 

with (9 W acting on the position s of the tensor product. The collective observables play the 
role of sufficient statistics for our model, and we would like to understand their asymptotic 
behaviour. Since all systems are independent and identically prepared, and the terms in each 
collective observable commute, we can apply classical Central Limit techniques to show that 
the following convergence in distribution holds for the collective states of" _ , 

qjjc(n) l 



Pjjcin) _£ 

yfn 
b,{n) £ 



tf(Re(z J y t ),Vj Ut ) 1 \<j<k<d; 
N(lm(zj,k),Vj ik ), l<j<k<d; 
N {u u Ai(\ - A,)) , \<l<d-l, 



yfn 

The key observation is that the unknown parameters M = (z, u) are recovered as means of the 
limit Gaussian distributions. However the limit model is not a classical one due to the fact 
that the collective variables do not commute with each other. Therefore we need to take into 
account the commutation relations of the limit variables by using the quantum Central Limit 



Theorem (CLT). These form a general continuous variable system as described in section 2.1 



with d 2 - 1 coordinates Z = (Z\, . . . , Z^_i) = (Qjj, P^k, B,) having the commutation relations 

[Z a ,Z b ] =Tr(p [O a ,O b T)l. 

By property (ii), this means that (Qj^,Pj,k) are canonical coordinates of d(d - l)/2 mutually 
commuting one-mode systems, and B, are classical random variables in the sense that they 
commute with each other and with the quantum coordinates (Qj t k, Pj,k)- By the CLT the limit 
state is the Gaussian 0(M, V) with mean M and covariance matrix V defined in property (iii), 
equation ( fT8| . In particular, the individual modes (Qj t k,P j,k, Bf) are independent of each other 
and of the classical variables B h and the latter have covariance matrix V cl , cf. < } 1 8] > . 

We define now the quantum-classical Gaussian statistical model f as 

9 := {<D(M, V) : M € R d ~ l x C J( ^ 1)/2 }. (20) 

and enunciate the Local Asymptotic Normality Theorem [19] which is used in establishing 
the optimality results in section[4] 

Theorem 3.1. The sequence of qudit models Q„ defined in ( |14| > converges to the quantum- 
classical Gaussian model V defined in \20\ , in the sense that there exist e > and channels 

T„ : M(C'f" -> 71 (L 2 (Rf d (flM)/2 ) ® L 1 ^ 1 ), 
S„ : T] (L 2 (R)® d(d - 1)/2 )®L 1 (R''" 1 ) -> M(C d f, 
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such that 



lim sup |p®» -S n (<p(M,\)) =0. (21) 
Note that the statement of the above theorem is different from that of the CLT in that Theorem 



3.1 shows that the collective state of the ensemble can be transferred by means of physical 
quantum channels to a Gaussian state, with vanishing norm-one error uniformly over the 
unknown parameters. The CLT is instead a statement about the convergence in law, for a 
fixed state and does not have an immediate operational interpretation. In the next section we 
will use LAN to transform the qudit reversal problem into a corresponding Gaussian reversal 
one, thus exploiting the solution to the latter obtained in section [2] 

4. Optimal channel reversal on mixed qudit ensembles 

We now focus on the main aim of this paper, that is to find optimal channels 

R„ : M(C d f" -> M(C d f m 

which reverse (at rate T = m/n) the action of the tensor product channel C m acting on n 
identically prepared qudits, cf. Q [see Fig.[T|a)]. The performance of channel reversal can 
be quantified by a figure of merit, or risk, given by the trace-norm error (^(p®") - p® m \ x . We 
adopt a frequentist approach and look to minimize the maximum risk over all input states. We 
actually work with a more refined version of the maximum risk known as the local maximum 
risk, which was already employed in other quantum statistical problems ifTTl [181 [341 l8l . For 
each state p, this is defined by 

Kmax(/? n ,P,r):= SU Pi \\Rn(l®")-T® m l 

l|T-p||<>ri +e 

and quantifies the worst performance of R„ in a n ^-neighbourhood of p. This restriction 
does not amount to making an assumption about the location of the unknown state, since 
using a small sample n l ~ e «nof the systems, one can effectively localize the unknown state 
within a confidence region of size n~5 +e . We aim to find an optimal reversing strategy whose 
asymptotic risk is equal to the local minimax risk 

7?minmax(P, O := lim SUp inf 9? ma x(^«,P, O . (22) 

The problem of finding the optimal reversing channel is thus reformulated as follows. For 
a given state p and output rate T = m/n, what is the minimax risk 9? m i nm ax(p> F) and which 
procedure achieves it? We will not solve this problem for all (p, F), but rather we will show 
that for each p there is an interval [0, T*] for which the reversal can be performed perfectly 
(i.e. with asymptotically vanishing risk 9? m inmax(p> O = 0), and we will subsequently identify 
an asymptotically optimal sequence of reversing channels R* which achieves the maximum 
rate F*. 

As anticipated, the quantum LAN theory IfTTl [T8l [T9l l20l |2D . reviewed in section [3] is the 
key tool in solving this problem. For technical reasons related to the validity of LAN, we 
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assume that the state p is not on the boundary of the state space, and is not degenerate so 
that its spectrum A - (A\,... ,Aj) satisfies A\ > ■ ■ ■ > A^. A typical random state satisfies 
these assumptions. Following equations (p"3j)-(fT6)) we can parametrize the states in the n~^ 2+e 
neighbourhood of p as p M / ^ A where M = (z, u) e ^ 2 x R d ~ 1 is a vector of expectations 
for the matrices \0\, . . . , 0^._\\. By LAN, the sequence of input quantum statistical models 



is asymptotically equivalent in the sense of Theorem 3.1 to the sequence of Gaussian models 



f„ over a quantum-classical continuous variable system with coordinates Z 

n:={0(M,V):||M||<« e } 



where <E>(M, V) which is a tensor product between d(d- 1)/2 independent one-mode displaced 
thermal equilibrium states $>(Zjjc, Vj,k) ( one f° r eacn P a ir j < k) with mean z^k and covariance 
vj^k (cf. equation (fl"8j)), and a classical (d - l)-dimensional Gaussian probability density 
JV(u,V cl ). 

We apply the same procedure to the output state p' := C(p) and denote by 

the transformed state belonging to the neighbourhood of C(p). The corresponding Gaussian 
model with coordinates Z' is 



P' n :={<D(M',V'): ||M'|| <n e ] 

where V depends solely on the spectrum of p' . The displacements of the two Gaussian models 
are related by a linear transformation 

X T : M i-> M' 



which describes the local action of C around p, and hence depends only on these two objects. 
The ensemble transformation C m acts on the localized states as 



C m : p 



(Ml ^ii,X) 



(M7 VV)' 



(23) 



Since both the input and the output and can be approximated by Gaussian states and the 
transformation between the displacements is linear, one can expect that the two Gaussian 
states are connected by a Gaussian channel G as depicted in diagra m The fo llowing 



theorem states that this is indeed the case and the proof can be found in Appendix A 
Theorem 4.1. There exists a Gaussian channel G such that 

G : <f>(M, V) i — > <f>(X T M, V). 

The matrices (Xc, Yc) characterizing G (cf. |5]lJ are Xg = X and Yc = V — X r VX. 

We now move to the problem of finding the optimal reversing channel for multiple qudits. 
As illustrated in the right side of the diagram (|3j, the solution relies on LAN to recast the 
problem into the Gaussian one solved in Theorem |3T| We now state the second main result 



of this paper, whose proof is in Appendix B 
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Theorem 4.2. Let p e M(C d ) be a full rank state with non-degenerate spectrum, and let 
C : M(C d ) — > M(C d ) be a quantum channel which is invertible as a linear map. Let (k*,H*) 



be the optimal rate and channel for the Gaussian reversal problem (cf. Theorem 2.1 



(f>(M, V) t - — > cf>(X T M, V) ^ cf>(kM, V) M e R d ~ l x C dW_1)/2 
where V, V, X depend on C and p as above. 

Then the maximum reversal rate for the qudit ensemble is F* = (k*) 2 and 

R* :=S r .„oH* oT„, 
is a sequence of optimal reversing channels i.e. 

limsup sup ||i?„(T®' I )-T® m [| 1 =0. 

l|r-p||i<fl 2 

This result shows, in full generality, that in order to reverse the action of an arbitrary channel 
acting on a large ensemble of n qudits, independent and identically prepared in an unknown 
mixed state p, the optimal strategy is to take a shortcut through phase space [35 1, optimally 
reverse the corresponding Gaussian channel, and then map the output back onto m qudits 
(with rate m/n = T*). 

This kind of digital-analogue interconversion might prove handy when it comes to counteract 
the effect of noisy channels affecting quantum states of, say, n-atom ensembles: The best 
strategy ought to be transferring their state onto light modes (e.g. using quantum non- 
demolition interactions |10|), implementing corrective procedures on the obtained Gaussian 



beams in phase space (according to the recipe of Theorem 2. 1 1, and then mapping the restored 
state back onto the atomic storage unit. We will now showcase our results on a specific 
example. 



5. Example: reversal of the phase flip qubit channel 

The qubit phase flip channel is defined by C(p) = E^pE^, with Eq = yfp Diag(l, 1) and 
E\ — s/1 - p Diag(l, -1) [ 1 1. Phase flip is a common source of error in quantum computation 



and several algorithms have been developed for its correction 1361 . According to Theorem 4.2 
an optimal way to reverse the composite action of individual phase flips on a n-qubit register 
is to map the state of the qubit register onto a Gaussian state (which in this case comprises 
a single quantum mode and one real classical random variable), reverse the corresponding 
Gaussian channel (up to a rescaling factor k*), and map the output back onto the qubits. 

As described in the previous section we consider states in a neighbourhood of p; the latter has 
in general the Bloch vector representation 

P=\{} + 1 rH (24) 



2 \ r x - ir y 1 - r z 

The neighbourhood of p is parametrized as in fP3j ) with M = (z, u) e C x R. The action of the 
channel C is 

C : r = (r x , r y , r z ) » r' = ((2p - l)r x , (2p - \)r y , r z ). 
From this we can compute the covariance matrices 

V = Diagf — , — ,1 -Hr|| 2 |, V = Diagf— , — ,1 
S \2||r||'2||r||' /' g l2||r'||' 2||r'||' 





Figure 3. (Color online) Plot of k* as a function of r x and r, for p = 0.6. 



where ||r'|| 2 = (2p - l) 2 (rj + r 2 ,) + r 2 . Similarly the matrix is obtained after computing the 
transformation of the local states under C 



X T - 

G - 



(2p- 



1) 















4(1- p)pr x r. []rj 
l|r'l|-||r|| ||r|| 



(25) 



We are interested in evaluating the threshold factor k* as given by (f8]l, which is associated to 
the optimal rate m/n = {k*) 2 for phase flip reversal on the qubit register. Due to the inherent 
symmetry of the channel, we may set r y = without any loss of generality, and assume 
I < p < 1. Therefore k* will depend on p, r x , r z . The results are summarized as follows. 

For p — 1 the channel is trivially the identity, therefore k* — 1 for any r. For p = 1 /2 the 
channel is never reversible, and k* = for any r. In the limit of pure input states of each qubit 
(|r| — > 1), we find that that k* — > 2p - 1 if r x — > 0, and k* — > otherwise (i.e., the channel 
is not reversible for any pure state, if there is a nonzero component on the x axis). In case 
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r x = 0, we find k*(p, r x — 0, r z ) — 2p - 1, which means that, as intuitively expected, the rate 
does not depend on r z . On the other hand, in case r z = 0, we have 

k*(p, r„ r z = 0) = ^[(2p - 1) 2 (1 - rl)]/[l - (-1 + Ipfrf]. 

We observe that for small values of r x this tends again to the line 2p - 1; for r x close to 1, k* 
is instead a sublinear function of p, as depicted in Fig. [2] showing that the reversal becomes 
inefficient in terms of the number of perfectly retrieved copies. In general, if a nonzero r x is 
present as well, the behavior of k* as a function of p is qualitatively similar to the previous 
case, but typically k* increases, at fixed ||r|| and p, with increasing ratio \rjr x \. The general 
shape of k* as a function of the Bloch vector components interpolates among the various 
limits discussed above, and is plotted in Fig. [3] 

6. Conclusions 

In this paper we have solved the general problem of optimal channel reversal for an ensemble 
of generally mixed independent and identically prepared qudits. This substantially extends 
our earlier work (8) on the particular instance of state purification (and dilution) of qubit 
ensembles. To accomplish our task, we have employed the versatile statistical tool of LAN 
IfTTl [19l to recast the problem in terms of Gaussian states and channels, then solved the open 
problem of optimal reversal of a general Gaussian channel acting on a quantum-classical 
Gaussian state with a given covariance. 

The methods demonstrated in this paper provide powerful strategies for counteracting 
undesired noise effects in quantum memories and long-distance quantum communications, 
based on interfaces between discrete ensembles and continuous modes |]9] [TUJ [T3] . A further 
generalization of our work, worth addressing in the future, could be that of considering 
the optimal reversal of quantum channels acting on correlated copies of n qudits, possibly 
requiring an extension of the LAN theory beyond the paradigm of independent and identically 
prepared systems. 
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Appendix A. Proof that the Gaussian Channel G in the diagram |3) exists 

Note that by the LAN construction the coordinates Z and Z' of the input and output Gaussian 
continuous variable systems have the same symplectic matrix Q., 

[Z a ,Z b ] = Tr(p[O a , O b ])l = iCl atb l = Tr(p[0' a , 0' b ]) = [Z' a , Z'„] (A. 1) 

and the covariance matrices are determined by the states p and p' = C(p) 



{Z a * Zb) = Tr(pO a * O b ), (Z' a * z;> = Trip' 0' a * 0' b ). 



(A.2) 
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Recall that any Gaussian channel from the input to the output is of the form 

G* : W, i — > Wx GZ e^ zTYaZ 

where Xc, Yg are real matrices of dimension d 2 — 1, with Yq positive and satisfying the matrix 
inequality 

Y G > i(X T G DX c - Q). (A.3) 

Since the means are transformed asMn X^M the matrix X G must be equal to X. Moreover 
the output variance is 

V = X r VX + Y c 

which determines the second matrix Yc = V - X r VX. It remains to show that Yc satisfies 
( |A.3| > which is equivalent to 

V + z'O > X r (V + /Q)X. (A.4) 
By using ( |A. 1 ] > and ( |A.2[ ) the latter can be translated into the following property of C 

V > X r VX (A.5) 

where \ a h = ti(C(p)0' a O' b ) and \ a j, = ti(pO a Ob). Finally, to verify this inequality note that 
for all M' 

Tr( P(IVM )X r 0) = X T C M = Tr(C(p(M^)0') = tx(p W )C(0')) 

which implies O' = X r O. This result together with the inequality T(A*A) > T(A*)T(A) valid 
for any channel T imply ( |A.5[ ): 

a,b a,b 

>£c^Tr(pC(0' fl )C(0p) 

a,b 

= c f X r VXc (A.6) 
Hence if C is a completely positive map then G is a Gaussian channel. □ 

Appendix B. Proof of Theorem 

We need to show that the reversal of C is possible on the interval [0, T*] and not possible for 
T > T* . Beginning with the qudit state 

rv n \ m — r, m 

L(P(M/ViU)) -P(M'i^,A>y 

where M' = X^M, we apply the LAN channel T n (cf. Theorem lOl followed by H* (cf. 



Theorem 2.1 1 to obtain a state which is asymptotically undistinguishable from the Gaussian 
state (5)(k*M, V). To this we apply the inverse LAN map S ,„ = 5V 2 « to achieve the corrected 
state 



7?maxOCp,r)= SUp 

[|M||«Sn« 
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We first show that R* is an (asymptotic) reversing channel for p m : 

ism _ d* I m \ 

^(M/vn,^) " V (M'/ vV)/ 

+ sup \\s m U(k*M,\))-S m oH* oT n (pf^ A J\ 

\\U\\<n' U ' 111 

< sup L(**M, V) - H* o T„ (p®" )|| + o(l) 
= ||*(**M, V) - ff* (0(M', V'))||, + o(l) = o(l), 



where we have used the contractivity of the trace norm under quantum operations, the 
properties of LAN and the fact that H* is a Gaussian reversing channel for G, cf. {7). 

We now show that it is impossible to obtain asymptotically exact reversal for F > F*. Indeed, 
suppose there exists a sequence of reversing channels R n : M(C d ) m — > M(C d )® Yn such that 

limsup sup ||^(p* )-p«*J =0. 

n->oa [|M[|<ra« 11 v " 111 

We then use LAN again to show that there exists a G-reversing channel T m o R„ o S„ with rate 
k > k*, which is impossible: 

\\T m oR n oS n ^(M',Y')-(P(kM,\)\l 

=0(1). 
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